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Abstract 

We study the sum (h(s) = ^Z^J S over the eigenvalues Ej of the Schr~dinger 

3 

equation in a spherical domain with Dirichlet walls, threaded by a line of magnetic 
flux. Rather than using Green's function techniques, we tackle the mathematically 
nontrivial problem of finding exact sum rules for the zeros of Bessel functions J„, which 
are extremely helpful when seeking numerical approximations to ground state energies. 
These results are particularly valuable if v is neither an integer nor half an odd one. 
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1 Introduction 



There is a large class of physically interesting problems where the Hamiltonian is closely 
related to the operator of Bessel's equation, this includes, e.g. : (i) classical equations 
for vibrating strings and drumheads, heat conduction in cylinders, normal modes in 
resonant cavities, Fraunhofer difraction through circular apertures, etc.; (ii) quantum 
particles which move freely within cylindrically or spherically symmetric domains, on 
whose boundaries the wave function vanishes (they play the role of reflecting or bouncing 
walls). The two-dimensional case has been made further complicated in two different 
ways: by altering the shape of the boundary — quantum billiards Jl|, |2| — and by intro- 
ducing magnetic potentials such that the domain is threaded by a line of magnetic flux 
— Aharonov-Bohm quantum billiards |}|. 

This paper deals with the extension of the last situation to spherical domains 
of arbitrary radius, a, diametrically threaded by a flux line. For such a system, the 
time-independent Schr" dinger equation takes the form 

J-[-iW - gA(r)] 2 M/(r) = E*(t), (1) 

where \i is the mass of the particle, q its charge and A(r) is some suitable vector 
potential, which produces a straight line of flux along the z-axis. To this end we choose 

A(r) = ^^e„ (2) 
2ir r sin 9 

where $ is the value of the flux and denotes the unitary tangent vector for increasing 
ip (the usual spherical coordinates, r,9 and ip, are employed). This corresponds to 
selecting a gauge in which the lines of the vector potential become concentric circles. 
After noticing that V • A = 0, as should be expected, ([T]) reduces to 

*(r,0 )¥ O = O, (3) 



r A sin 9 \ oip , 



where k 2 = ^ 2 and a = ^ — . Applying the usual variable separation method, we set 



r h 2 

ty(r, 9, ip) = R(r)y(9, <p) and arrive at 

R{r) = 0, (4) 



|-r 2 ^ + Pr 2 - C 
or or 



2 



1 _ 2m d_ _ A + c 

sin 9 89 " 39 ' sin 2 9 \ dip 2 dip J 



i a a 

sin 9— 



0. 



(5) 



with C a separation constant. Writting C = (3{f3 + 1), (fjj) becomes the spherical Bessel 
equation of parameter /3, not necessarily integer. Thus, the solution regular at r = is 



R(r) ~3f3(kr) = J^Jp +1/2 (kr) : 



(6) 



where jp denotes the spherical Bessel function. As a result of the boundary condition 
imposed on the surface r = a, ka must be a zero of J/3+1/2 ■ The energy is therefore 
discretized in the way 



E, 



n 1 



8n 



2fia : 



~3f3+l/2,ni 



(7) 



jf3+i/2,n being the nth positive zero of Jp+i/2- 

Taking y(9, ip) = Q{9)<f>((p) with the function 4>(<p) of the form e imip , m G Z, and 
doing the change of variables x = cos 9, eq. (El) yields 



d 



d 



l-x 2 )—-2x— + (3((3 + l 



(m — a)' 



e(9(x)) = 0, 



dx 2 dx 1 — x 2 

and writting Q(9(x)) = (1 — x 2 y m - a ^ 2 u(x), the function u(x) satisfies the following 
equation 



(1 - x 2 ) u" - 2x (\m - a\ + 1) v! + [C - \m - a\(\m -a\ + l)]u= 0. 



(9) 



Since x = is a regular point, we make the assay u(x) = ^2 a kX k , which leads to the 

k=0 

recurrence 



(k + \m — a\)(k + \m — a\ + 1) — C 
a k+2 = — — , — — °fc- 



(10) 



(k + 2)(k + 1) 

Good behaviour at x — ±1 can only be achieved by truncating the series into a polyno- 
mial. This takes place if there is some p G N such that C = (\m — ot\+p)(\m — at\+p+l). 
This recurrence would be the same as the one for Legendre's associated functions P/", 
were it not for the presence of a (generally) noninteger a. Therefore, we adopt the 
notation 



6(0) ~ P^- a (cos9), \m-a\< p, 



11) 



where (3 = \m — a\ + p is also non-integ rule. The general form of the wave 

function is a superposition of the sort 



n,Lm 



*(r, M = £ J/3 J/3+iAn- P^costf) e*"* (12) 



The angular eigenfunctions ~ P™~ a (cos6) e im<p generalize the spherical harmonics 
Yi m (6,(p), corresponding to a = 0. It is precisely for this reason — i.e. the presence 
of a — that the solutions are no longer eigenfunctions of the angular momentum L 2 
(they are instead eigenfuctions of the mechanical momentum II 2 ). This was to be ex- 
pected, as radial symmetry has been broken by A(r). The energy spectrum has now 
turned into 

h 2 

Ep, n = 2^2^?+i A»' m G Z ' P = \ m - a \ + Pi V e N. (13) 

The ground state energy is therefore = Y^jf a \+i/2,v ^ or \ a \ — V^- Notice that 
a = 1/2 gives rise to double degeneracy, since states with m and — m + 1 have the same 
energy. This is also true, in particular, for the ground state. 

Let us recall that the zeta function of a Hermitean operator A with infinitely 
many discrete eigenvalues Aj is defined as 

Us) = £ ^ = TrA- (14) 

An alternative to the computer evaluation of these Bessel function zeros is their 
approximation by means of sum rules. When A has a discrete spectrum Ai < A2 < 
A3 < . . ., several ways to find estimates of low-lying eigenvalues — quoted in @] — are 
in order. The most adequate here is Euler's method (originally applied to the zeros of 
Jo), which is based on the inequalities 

Us)- 1/S < A, < -4^tt- (15) 
(a{s + 1) 

C (^) 

In fact, lim Ka(s)1 _1 ^ s = lim ——, r = Ai and, prior to taking the limit, (PTSI) always 

s ^°° s ^°° ( A (s + 1) ' ' 1 — v 

holds. 
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In our case, we will take A to be the operator of the spherical Bessel equation 
(§J) with C = (\m — a\+p)(\m — a\+p+l). Thus, for each possible I, the zeta function 
Ca{s) to be considered is of the type Ca(s) ~ ( u (2s), with 

OO 1 

Us) = E — » Re s > ( 16 ) 

n=l ^^n 

where j vn is the nth nonvanishing zero of The condition Re s > 1 stems from 
considering the asymptotic form of the j vra 's, which is roughly (v — ~J ~ + 7m, and from 
comparing with the Hurwitz zeta function. At any rate, the existence of poles at some 
values of s ought to be checked by operator heat-kernel expansion methods. 



2 Recursive Rules for Zeta Functions associated to 
Bessel Operators 

2.1 Quadratic Recursion 

When written in the form of an infinite product involving its nonvanishing zeros, the 
z/th Bessel function of the first species reads 

z v 00 / z 2 \ 

Differentiating the log of this expression with respect to z, and using the identity 

J'Az) = -J v {z)-J v+l (z), (18) 

z 

we arrive at 

W = (19) 

1 °° f z \ 2m 

and, since j v \ < j u2 < taking \z\ < j v \ we are able to put k = — ) 

\3mJ 

Then, interchanging the summations and looking at ([16]) we come to 

%^ = 2 fc,(2n)^- 1 . (20) 



This formula will constitute our starting point and is a well known result (see e.g. 
|§, pg. 61). Now one can, of course, find the values of the zeta functions through 



C(2n) 



1 d 2n - x J u+1 {z) 



2(2n-l)! dz 2 ^ 1 J v {z) 



2=0 



. Rather than doing this, we take the derivative 

u+1 



of (PH) with respect to z, and making use of ( |18D and of J' u+1 (z) = J v (z) '■ — Ju+i( z ), 



we obtain 



JJz) 



JJz) 



£C/(2n) 2(2n-l) 



„2n-2 



(21) 



n=l 



Replacing Jv +^ with the r.h.s. of fl20|), we get a sum of series which can be recast into 
the form 



i + E 



71=1 



71-1 



-4(i/ + 7i)£,(2n) +4^ C.(2/)C(2n - 2/) 



i=l 



^ 2 "~ 2 = 0. 



(22) 



Since it must be satisfied for any z with \z\ < j u i, the identity must hold separately for 
every power of this variable. In particular, for n = 1 we are led to 

1 



C(2) 



4fi/ + 1 



while, for n > 2, the 'quadratic' recursive relation 

1 71-1 



G(2n) 



n + v 



X;C,(20C,(2n-2Z) 



(23) 



(24) 



2=1 



follows. In Appendix A we present an alternative derivation of this formula by direct 
integration over a convenient contour in the complex plane. 

This equation, together with the previous, provides the value of C, v {2n) at any 
even integer argument. The first results thus obtained are: 

1 



C(4) 
C,(6) 
C*(8) 

C,(io) 



2> + l) 2 (z/ + 2)' 
1 

2 5 (z/+l) 3 (z/ + 2)(z/ + 3)' 

5^ + 11 

28(i/ + l) 4 (z/ + 2)V + 3)(z/ + 4)' 
7z/ + 19 

2> + \f[v + 2) 2 (i/ + 3)(i/ + 4)(i/ + 5) ' 



(25) 
(26) 
(27) 
(28) 
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These expressions are quoted in || as due to Rayleigh. Now, with the help of (|24|), their 
derivation has been straightforward. 

An immediate consequence of (p4|) is the nontrivial equality 



n— 1 oo i 

E E -W-^I = (- + l)C„(2n). (29) 

/—I fc,m=i JvkJvm 

The procedure here described finds application to other functions that can be 
expressed as infinite products. For instance, we take sin z = z TT ( 1 — — . Its 



n=l 



^ 2 



logarithmic derivative gives (assuming |z| < n) a known expansion for ctg z: 

i oo ^2m+l 

ctg 2 = --2E((2m + 2)— ^, (30) 

^ m=0 " 

where ( is the Riemann zeta function itself. Differentiating with respect to z, the l.h.s. 
will be — esc 2 z = —(1 + ctg 2 z). Replacing ctg z with the series (p0[), and the r.h.s. 
with its derivative, we obtain an identity which must be fulfilled for every power of z. 
Thus 

' C(2) 6 ' 

1 n-l (31) 

C(2n) = — -iX)C(2n-20C(20, ™ > 2- 

n + 2 i=l 

This provides a method for the recurrent evaluation of C(2w) from just the knowledge 

of C(2)- Bearing in mind the relation between zeta function and Bernoulli numbers, 
(27r) 2n 

C(2n) = (~l) n ^, -B2n; the previous equations read 

h 

(32) 

£ 2n = "TT— T E I I B 2n-2iB 2l , n>2. 

They supply a different way of finding all the nonvanishing Bernoulli numbers {B\ aside) 
from the value of B 2 . Eq. fl3ll) yields also 

Jl— 1 OO -r i q 

S£^^ CW ' (33) 

Actually, all these results can as well be obtained from those for by taking v — |, 
since Ji/ 2 (z) = W^sinz. 
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2.2 Linear Recursion 

The above recurrence has turned out to be a powerful tool for successively obtaining the 
zeta functions of even argument. But, having succeeded in finding a 'second order' law 
— as the r.h.s involves products of zeta functions at lower arguments — we would also 
like to have a linear rule available. The example of Euler's method (shown, e.g. , in JIJ, 
pg. 500) takes advantage of a lin ear recurrence among the £o(2n)'s. Now, our aim is to 
find the general form of such type of relation, thus extending the procedure to any v. 
We start, once more, by taking fl20"D , this time written as 

oo 

J u+1 (z) = 2J v {z) £ U2n + 2)z 2n+1 . (34) 

n=0 

Next, replacing the J v 's with their power series expansions 

*W = £ TT& fa) (35) 



k\(k + v)\ \2, 

(in general, these factorials are to be understood as gamma functions (k + v)\ = T(k + 
v+ 1)), we obtain a relationship between series which, after some index rearrangements, 
reads 



E 



n=0 ^ 



feo2fe 



1 n (-l) fe 2 

2 n!(n + i/+ 1)! ~ 2 £ (n - fc)!(n - fc + z/)!^ 2fc + ^ 



z2n+v+l = Q ^ 



Validity for any z calls for the vanishing of every coefficient. Therefore 

i (^ TT y - g'-^ V'qUn + t-V ' 2 * + 2) ' " " °' (3?) 

which amounts to 

1 ( n — l\ ( v + n — 1 \ 

£(-l) fc 4W C.(2fc + 2),n>l. (38) 



4(r/ + n) ^ \ jfe J \ k 

This is the linear relation we were looking for. 
The first resulting identities are: 



1 



4(i/ + 2) 
1 

4(i/ + 3) 
1 

4(i/ + 4) 
1 

4(i/ + 5) 



C„(2)-4(i/ + l)C„(4), (40) 
£,(2) - 8(1/ + 2)&,(4) + 32(i/ + l)(z/ + 2)C„(6), (41) 
C„(2) - 12(i/ + 3)C(4) + 96(i/ + 2)(i/ + 3)C„(6) 

-384(z/ + 3)(z/ + 2)(z/+l)C(8), (42) 
C„(2) - 16(i/ + 4)C„(4) + 192(i/ + 4)(i/ + 3)C(6) 
-1536(i/ + 4)(i/ + 3)(u + 2)C„(8) + 6144(i/ + 4)(z/ + 3)(z/ + 2)(z/ + 1)C„(10). 

(43) 

This can be regarded as a system of linear equations in C/(2), Ci/(4), . . . , C(10), which can 
be solved by repeated substitution from the first. It can be checked that its solution is 
(p3|) , fl25|), (p6|), (|27|) , (p8|) . The following remark is in order: although a linear recursion 
looks in principle simpler than one involving products, in practice it is much easier to 
work with (1241) than with (B8T) . 



3 Numerical evaluation of the ground state energies 



Euler's method, combined with (|23"D , (|24]) , turns out to be a very efficient procedure for 
the computation of j v> i, associated to the ground state (here v = \m — a\ + p + 1/2). 
The lower and upper bounds in (|15D are now [( u (2k)]^ 1 ^ 2 ^ and ( u (2k) / ( u (2k + 2), 
respectively, where the square root comes from the factor 2 in the argument when going 
from to (a- We calculate them for increasing fc's until convergence of both the lower 
and the upper successions becomes apparent. At every step, a new ( v (2k) is found by 
means of a recursive function — implementing ( p3[) and (p4 ) — which is part of a simple 



C program. Table |T| shows the intermediate results obtained for the cases v = 0, 1/2 
and 1, whose first nonvanishing zeros are known to be 2.404826 . . ., 7r and 3.831706 . . ., 
respectively. The figures obtained after arriving at k = 10 for different z/'s between 
and 1 are displayed in the right box. Of these results, only those for v > 1/2 can be 
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physically meaningful, as the lowest possible ground state corresponds to v — 1/2. Such 
value would be attained only for a = 0, while higher z/'s correspond to nonvanishing 

a's. 

The calculation has been repeated — with the same accuracy — using a program 
which implements the linear recursion instead of the quadratic one. No figure of the 
results shown in Table 1 has changed, alhough the execution time is now shorter. How- 
ever, away from this range the numerical errors produced by both algorithms may differ. 
Actually, for v close to 100, we have observed differences of the orders of 10~ 6 and even 
10~ 5 within the first ten steps. 



k 


V 


= 


v = 


1/2 


V 


= 1 


1 


2.000000 


2.828427 


2.449490 


3.872983 


2.828427 


4.898979 


2 


2.378414 


2.449490 


3.080070 


3.240370 


3.722419 


4.000000 


3 


2.401874 


2.412091 


3.132603 


3.162278 


3.812737 


3.872983 


4 


2.404424 


2.406133 


3.139995 


3.146427 


3.827710 


3.843076 


5 


2.404766 


2.405069 


3.141280 


3.142768 


3.830778 


3.834980 


6 


2.404816 


2.404871 


3.141528 


3.141883 


3.831478 


3.832667 


7 


2.404824 


2.404834 


3.141579 


3.141665 


3.831648 


3.831990 


8 


2.404825 


2.404827 


3.141590 


3.141611 


3.831691 


3.831791 


9 


2.404826 


2.404826 


3.141592 


3.141597 


3.831702 


3.831731 


10 


2.404826 


2.404826 


3.141593 


3.141594 


3.831705 


3.831713 



V 


M 





2.404826 


0.1 


2.557451 


0.2 


2.707073 


0.3 


2.854097 


0.4 


2.998849 


0.5 


3.141593 


0.6 


3.282545 


0.7 


3.421890 


0.8 


3.559780 


0.9 


3.696347 


1 


3.831705 



Table 1: Left: successive values in the approximation to corresponding to v = 0,1/2,1. 
For each is, the values on the left and right columns are the lower and upper bounds, obtained 

as [Cu(2k)]' 1/(2k) and [( u (2k)/( u (2k + 2)] 1/2 respectively. Right: different >i's for some values 
of v between and 1. 
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4 Conclusions 



We have focused in this paper on the derivation of exact and explicit formulas for the zeta 
function corresponding to the zeros of the Bessel function, ( u (2k) (for arbitrary u), and 
in their systematic use for obtaining estimates of ground states of related Hamiltonian 
operators. The accuracy of the procedure here developed has proven to be excellent. 
Yet, the main advantage of the method has not been exploited fully. By this we mean 
the possibility of finding such approximations even in cases where the spectrum does not 
emerge from standard special functions, but can nevertheless be calculated, at least 
for some special values 0. 

Another subject to be explored is the deformation of the spherical domain. The 
corresponding case in two dimensions was tackled by means of conformal mappings be- 
tween the unit disc and other noncircular domains @ . In three dimensions the difficulty 
would surely increase because such powerful complex transformations are lacking. We 
think that these subjects are worth studying in further detail. 

In particular, here, by dealing with a problem of known spectrum, we have been 
able to relate the physical information concerning the energy levels with some remarkable 
numerical properties involving Bernoulli numbers and ordinary Riemann zeta functions. 
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A Appendix: Derivation of the quadratic recurrence 
by complex integration 



The zeta function 

oo 

can be written as the integral 



m=l 



CO) 



/c*xi^/c* ,n i J «i* 



-2-1 



(44) 



(45) 



2ni Jc J v (t) 

where C is the sector-like contour of the complex plane domain |arg(£)| < 9, with 
< 9 < ir/2, and \t\ > e, with < e < Ai and Re t > 1. Recalling that, for 
| arg(w)| < 7r — 5 (5 > 0) and \w\ — >• oo, 



J v (w) 



TCW 



Pyiw) COS ( W — (h> + 



Q u (w) sin — (z/ 



2 ; 2 



(46) 



where P v and are series of negative powers in w (we assume that v > —1/2), and 
particularizing it for the case w = pe ±l7T ^ 2 (p > 0), after substitution in expression (|45|), 
we obtain 



iZ 



2tt Y£Jc ±w 



/ 2ntvX Y ,(± / (t-^(u+^)) ) ,/(/) 



+ 2 



le 



-2 + 1 



1\ 7T 



(47) 
(48) 



l-z V 2/ 2z 

where we have taken here 9 = it/ 2 and C±(e) are the two halves of the small circular 



part of the contour, of radius e. The two integrals in ( |48D are analytical functions of 
z, while the last term is an explicitly meromorphic function. We thus see that ( u has 
simple pole at z — 1 with residue 

Let us now consider the restriction of to the strip — 1 < Re z < 0. The limit 
e — > yields the fundamental expression 



a,) = -sm( v -,j j / o rfpp- 



2-1 



111 



27rpe- p /,(p) 



-1 < Re^ < 0. 



(49) 
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It can be continued — in the value of Re z — to the right (and also to the left) of the 
real axis. Adding and subtracting in the integrand a term like e~ p J2t=o P k ( a k + bk In p), 
with ao = ln[\/27r/(2 I V!)], bo = z/+ 1/2, and adjusting the remaining a's and o's so that 
the difference with the integrand in be 



In 



2%pe- p I v (p) 



2p 

e-f J2 P*K + b k \ap)~0 (p 2p+1 lap) (50) 

k=0 



when p — > 0, it can be proven that the function G v (p) = In [2 v u\p u e p I u (p)} satisfies 

G" + G' 2 + G" ^2 + + ~~~p~ = °' ( 51 ) 

which has a solution of the form G(p) = J2T=i c kP k , c i = — 1; i n a neighborhood of 
p = 0. In terms of G u (p) and H u (p) = e p G u (p), the condition (|50D is written as 

2p 



H»(p) - £ a fc p fe + (e> - 1) In + M„(p) hip ^ (p 2p+1 lap 



(52) 



where M u (p) is analytical around p = and contains all the information on the b 
coefficients — which are rigorously proven to exist but will not interest us here. The a's 
are read off from 05lq), with the result 

k 



2ty 



Q,k = — In 

k\ \ 2 v v\ 



+ E 



(53) 



Back to (|49|), we see that 



— sm — z / ap p 

7T V2 yio FF 



2p 



2-1 



111 



27rpe- p I u {p) 



2p 

e~ P J2( a k + bk In p)p h 

k=0 



+ -sia(^z)Y / [^{k-z) + b k r'(k-z)} 



(54) 
(55) 



with the a's and 6's determined as above, provides the desired analytical extension in 
the strip — 1 < Re z < 2p + 1. After some straigthforward manipulations we see that 
for z a positive even integer, we get 



C(2m) 



l) m+ W 



(56) 
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It is also clear that the 6's do not contribute to the expression in the limit and, by 
working out the precise value of the c's from (|52|), we obtain C„(2) = [4(z/ + and 
the quadratic recurrence 

C„(2(p + 1)) = — ^— r E U2(* + l))Cv(2(p - k)), (57) 
valid for any integer p, p > 1. 
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